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A system of identical bosons with short-range (contact) interactions is studied. Their motion is 
confined to one dimension by a tight lateral trapping potential and, additionally, subject to a weak 
harmonic confinement in the longitudinal direction. Finite delay time associated with penetration 
t-H | of quantum particles through each other in the course of a pairwise one-dimensional collision in 

the presence of the longitudinal potential makes the system non-integrable and, hence, provides a 
mechanism for relaxation to thermal equilibrium. To analyse this effect quantitatively in the limit 
of a non-degenerate gas, we develop a system of kinetic equations and solve it for small-amplitude 
q , monopole oscillations of the gas. The obtained damping rate is long enough to be neglected in 

y-^ • a realistic cold-atom experiment, and therefore longitudinal trapping does not hinder integrable 

' dynamics of atomic gases in the ID regime. 

(N : 

^ I. INTRODUCTION 

Ctf ■ 

One-dimensional systems (ID) [l|, HJ often provide suitable models to study the fundamental processes of quantum 

■ dynamics, coherence and noise in interacting many-body systems. Integrability is one of such fundamental questions 
most closely related to the ID physics. In an integrable system [1,0], the number of integrals of motion is exactly equal 
to the number of degrees of freedom. Thus the observables, which can be measured experimentally, namely, single- 

q-i particle (momentum and position distribution) or few-particle (correlation functions) ones retain the information on 
the initial conditions forever. As a result, a quenched (i.e. suddenly driven off the equilibrium) integrable system 
may undergo relaxation towards the generalized Gibbs (or fully constrained thermodynamics) ensemble @, instead 
of the thermal equilibrium state. A similar type of relaxation may be observed even in non-integrable systems [f|. 
Also local quantities of the Bose-Hubbard model are shown [7| to approach the thermal Gibbs state much faster 
than non-local ones, such as the density-density correlation function. For a very wide class of non-integrable systems, 
however, the eigenstate thermalization hypothesis Q (see also @) may hold, enabling fast relaxation to thermal 
equilibrium. Of corse, the term "thermalization" with regard to closed (isolated from environment) systems means 
that all experimentally measurable observables approach their thermal distributions; the information on the initial 
conditions becomes contained in many-particle correlations (ultimately, in correlations involving all particles in the 
system), whose measurement is not feasible in practice. 

■ An uniform ID system of indistinguishable bosons interacting with each other via pairwise delta-functional po- 
! tential is known to be integrable and described by the Lieb-Liniger model [lol fill ]. It can be, up to certain extent, 

t^- ■ implemented in an ultracold-atom experiment in optical lattices [l2|, EH or on atom chips [T3 - [l7{ , provided that both 
the temperature and mean interaction energy per atom are well below the excitation energy of the first excited state 
. of the radial trapping Hamiltonian. Under these conditions, the radial motion of atoms is confined to the ground state 
7-H ' of the radial trapping Hamiltonian. Ultracold atomic systems deeply in the ID regime can be prepared in optical 
lattices [l8[ and on atom chips [l!|. However, in reality no system is perfectly ID, but the actual question is, on 
which timescale it can be described as ID. For ultracold atoms under tight lateral confinement one-dimensionality 
and, hence, integrability are lifted by atomic interactions causing virtual population of excited radial modes, even 
if the energy of thermal motion is too low, and thermal excitation of the radial modes is strongly suppressed by a 
Boltzmannian exponential factor. The role of the virtual radial excitations in the dynamics of ultracold atomic gases 
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in tight w aveg uides has been first noticed in the context of macroscopic flow of degenerate atomic gas through a 
waveguide [2(| and decay of mean-field solitons [2l|. Microscopically, virtual radial excitations give rise to effective 
3-body elastic collisions, which have been suggested as the source of thermalization in ultracold atomic gases on atom 
chips |22| (see also Refs. [UHil)- Interestingly, if the interatomic repulsion is so strong that the system approaches 
the Tonks-Girardeau limit, interatomic correlations reduce the rate of thermalizing collisions [23l425l |. The latter 
conclusion is in agreement with the results of the quantum "Newton craddle" experiment [l8j . 

However, this is not the only source of non-integr ability. An ultacold atomc cloud in a real experiment is also 
confined longitudinally, although the ratio of the fundamental frequencies of the longitudinal (oj||) and transversal 
trapping may be as small as uj^/uj± ~ 10~ 3 . Obviously, an ideal gas in ID in the presence of harmonic confinement 
is an integrable system, and so is the gas of impenetrable bosons [261 ] . However, the longitudinal confinement should 
lift the integrability of the Lieb-Liniger model for any finite atomic interaction strength. Estimating the corresponding 
relaxation rate is the subject of the present paper. 
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II. PHYSICAL MODEL 



The Hamiltonian of the system of N identical bosons of mass m confined to ID and also harmonically trapped in 
the longitudinal direction reads as 



N 



r^&_ mojj 2 h 2 c N 

2m dx 2 2 3 m 

3 i=i+i 



(1) 



Here c is the pairwise interaction strength. In the limit — > we recover the Hamiltonian of the Lieb-Liniger 
model El • Virtual excitations of the radial modes [20l - |2"4j are not included into the model Eq. ([1} , therefore 
we can investigate the integrability breakdown solely due to the longitudinal confinement. The latter is assumed 
to be harmonic in order to distinguish the effect under study from dephasing of the atomic motion occuring in an 
anharmonic trap. 

In what follows we consider a non-degenerate, weakly interacting ID bosonic gas. The non-degeneracy means that 
the atomic motion can be viewed as quasi-classical propagation in the longitudinal direction of distinct wave packets 
localized of a size much smaller than the mean interatomic distance. However, as we noted in Sec. HI the radial motion 
of atoms remains quantized, confined to the ground state of the radial trapping Hamiltonian. This is exactly the 
ponit, where our models differs from previous calculations of damping of excitations in a classical atomic gas in very 
elongated, but still three-dimensional (3D) traps 27]. In the latter case, (i) time delay in the process of scattering is 
neglected and (ii) collisions are essentially 3D and lead to scattering to a whole manifold of radial modes. 

A collision of two atoms is not an instantaneous event, but entails a time delay in propagation of the two colliding 
wave packets. The delay time is [28| 

u dk r ' ^ 

where u is the relative velocity of the collision (we assume here wave packet propagation in ID with u > 0), hk r = mu/2 
is the corresponding momentum (m/2 is the reduced mass) and (p is the phase shift of the transmitted wave after ID 
scattering on the given pairwise interaction potential. Having the formula tp = — arctan[c/(2fc r )] for scattering on the 
contact (delta- functional) potential [10], it is easy to obtain 

t = mc/[hk r {4:k 2 + c 2 )]. (3) 

Since in a weakly-interacting gas k r 3> c for almost every collision, 

r « mc/{Ahkl) = 2h 2 c/{m 2 u 3 ). (4) 

Formally, Eq. (jj) coincides with the classical expression r = [2/(mu 3 )} /_ °° dx Ui2(x) for the delay time of a fast 

particle of mass m/2 traversing the potential Ui2(x) ["fast" meaning that the particle kinetic energy mw 2 /4 at x — > ±oo 
is much larger than the maximum value of U\2{x)\. Of course, in a case of the repulsive (c > 0) delta- functional 
potential the classical description fails, but Eq. (j4]) holds also in the regime of quantum tunneling through a weak 
potential barrier. If the atomic interaction is attractive, we formally obtain r < 0. However, the theoretical approach 
developed in Sec. |TTT] assumes r > 0. Therefore we consider in this paper only gases with repulsive interactions (c > 0). 

The finite delay time docs not change the statistical distribution of atomic velocities in an infinite waveguide in 
the thermodynamic limit. However, the situation becomes different for the trapped gases. The relative motion of 
two particles is assumed to be classical in the intervals between collisions, and the center-of-mass motion is always 
treated classically. Assume that at time to two particles with the velocities t>i(io) = v c \ and Vzita) = v C 2 collide 
at the point X\(t) = X2(t) = x c . We assume the following model: during the collision time r, which depends on 
the absolute value of the relative velocity \v cl — w c2 |, the two atoms form a composite particle of mass 2m moving 
in the potential muj 2 X 2 , X being the coordinate of the composite particle, with the initial conditions X(to) = x c , 

V(to) = X(to) = (v c i +v C 2)/2. The kinetic energy of the relative motion is ascribed to an "internal degree of freedom" 
of the compound particle, thus ensuring energy conservation. After the time r the compaund particle disintegrates, 
releasing two atoms at the point 

u , ^ v c i + v c2 . 
£1,2(10 + t) = x c coscj||T -I smw||T. (5) 

II 

The velocities of the particles are determined by the center-of-mass velocity and the relative-motion kinetic energy 
stored internally in the compound particle and released back at t = t + r: 

,. , s ■ v cl + v c2 v c i - v c2 
vi,2(to + t) = —ojx c smcj||T H coso>||T ± . (6) 
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Since the atoms are identical, there is no one-to-one correspondence between the upper (lower) sign in front of the 
last term and the index (1 or 2) numbering the atoms. Although the total energy of the pair of the atoms is conserved 
exactly, the single-atom energies are not. They will be shifted by Ty( Vc1 ~ v c2)[\{v c i+v C 2)(l — cosu\\t)+l)\\x c sin w\\t]. 



III. KINETIC EQUATIONS 



To quantify the single-atom energy redistribution noticed in the previous section, we derive a set of kinetic equations 
describing ID atomic collisions. We assume a two-component model: the system consists of single atoms chrarcterized 
by a joint co-ordinate and velocity distribution function f(x, v, t) and diatomic compound particles representing atomic 
pairs, not actually bound but remaining close to each other during the delay time r, which is considered as the effective 
collision duration. The distribution function of compound pparticles F{x, V, u, t) has, apart from the time t, three 
arguments: x stands here for the compound particle co-ordinate (center of mass of two colliding atoms), V is the 
center-of-mass velocity, and u is the relative velocity of two atoms before the collision; u can be considered also as a 
label for "internal excitations" of a compound particle. Because of idistinguishability of the atoms, it is sufficient to 
consider only positive values of u. We assume that the gas is dilute, i.e., only a small fraction of the atomic ensemble 
is simultaneously contained in compound particles. A rate of ID collisions experienced by an atom with the velocity 
v is 

Ho dv ' 

\v — v'\ fix, v, t)j{x, v', t), and each collision is followed by compound particle formation (it is impossible 
for two identical, interacting bosons to tunnel through each other in ID without delay in propagation). 

Also we assume that a compound particle does not experience collisions with other atoms during its lifetime. To 
make the set of kinetic equations simple, we assume (instead of a well-defined (deterministic) time interval between 
formation and disintegration of a compound particle) that a compound particle undergoes exponential decay with the 
rate l/r„. The lower index index u explicitly indicates the dependence © of the delay time on the relative velocity u 
of the atoms before collisions. The products of the decay are two atoms emerging at the same point with the velocities 
V ± u/2. The equality of the relative velocities before and after the collision is specific for motion in a harmonic 
potential, where the center-of-mass and relative-motion degrees of freedom are decoupled. 

The corresponding set of kinetic equations now can be easily written: 

d d d 

— f{x, v, t) + v—f(x, v, t) - w\\x—f(x, v, t) = 



dv' 



\v - v'\ f(x,v,t)f(x,v',t)- 



^-F(x, V -±^,\v-v'\,t 

\v—v'\ \ * 



(7) 



—F(x,V,u,t) + V— F(x,V,u,t)- 







1 



w a— V, u, t) = Fix, V, u, t) + 

11 oV T u 

uf(x,V+^t)f{x,V-^,t). 



(8) 



It belongs to a broad family of equation sets describing chemical reactions on the kinetic level [2£ 
total number of atoms 



It conserves the 



N = 



dx 



dv 



f(x,v,t) 



du Fix, v, u) 



(9) 



and the total energy of the system 



E 



dx 



dx 



dve a f(x,v,t) + 

i 

dV / due cp F(x 1 V 1 u^t), 
Jo 



(10) 



where e a — miujf,x + v )/2 and s cp = m(wnj; + V~ + u /4) are the single-atom and compound-particle energies, 



respectively. Eqs. (jTl [Sj) are also consistent with Newton's second law d 2 Xj<sjdt 2 
Atv = dx x dv [fix, v, t) + 2 J °° du Fix, v, u, t)]/N of the iV-atom system. 



w; Xn for the center of mass 
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Eliminating F(x, V, u, t), we arrive, by an exact transformation, to an equation for the distribution functions of 
atoms only: 

d d d 

— /(a;, v, t) + Vg^f( x , v , t) ~ u \\ x g^f( x ^ w > f ) = 



dv' \v — v'\ 

oo 

g-«/TI„_„/l 

ds f(x s ,vj,t - s)f(x s ,v B ,t-s) 

T \v—v'\ 



f{x,v,t)f(x,v',t)- 

(11) 



where x s and vf 1 are given by the right-hand side of Eq. ([5]) and Eq. ([6]), respectively, with the following substitutions: 
x c = x, v c i = v, v C 2 = u', and r = — s. If we approximate the exponential kernel in the right-hand-side of Eq. (| 1 1|) 
by S(s — 77), 77—7- 0+, this right-hand-side (the ID collisional integral) becomes zero, and we obtain in this limit an 
equation of motion of a collisionless gas. The Maxwell-Boltzmann distribution 

Nmuju ( munx 2 mv 2 \ 
/o(*,«,t) = ~^jT exp y 2T~ J (12) 

is a stationary solution of Eq. (fTTj) . Here we set Boltzmann's constant to 1 and thus measure temperature T in energy 
units. We normalize /o in Eq. (|12p to the total number of atoms AT, since the interactions are week and the mean 
number of compound particles simultaneously present in the system is much less than N. 

IV. MONOPOLE MODE 

If a system state deviates from equilibrium only slightly, the typical time of return to the equilibrium is estimated 
from the linearized form of the kinetic equation [30j . Following the standard linearization scheme, we represent 
f(x,v,t) = fo(x,v,t) + fx(x,v,t) where the small perturbation term fi describes the deviation from the equilibrium 
and is orthogonal to /ojT^ dx dv fo(x, v, t) fi(x, v, t) = 0. For the sake of defmiteness, we consider one specific 
(monopole) ID mode [3l|, and choose the perturbation term in the form 

fi(x,v,t) = f (x,v)[A(t)(v 2 - u 2 x 2 ) + B(t)u\\xv]. (13) 

If we completely neglect collisions in Eq. (jlip . we obtain the following linearized equations: 

A + 2w||B = 0, S-2w||A = 0. (14) 

In the next iteration, we substitute the solutions of Eq. (fT4")l A(t — s) = A(t) cos2w||S + B{t) sin2w||S and B(t — s) = 
— A(t) sin2w||S + B(t) cos2w||S into the right-hand-side of Eq. (fTTT) and finally obtain 

A + (2w|| - (3)B + 7A = 0, B-2w||A = 0, (15) 

where 

Nunm? f°° . 2,,.^ 2uj\\t u 

(3 = -—Kpr / duu 5 e- mu -f-^ 16 

64ttT 3 J q 1 + 4w2 t 2 



and 

7 ~~ 64ttT 3 J """ " ~ 1 + 4cj|t, 

Solution of Eqs. (|17[) yields a complex eigenfrequency. Its real and imaginary parts give the oscillation frequency 



^ r d uu'e-™''W 4 "" T " (17) 



u = 2w|| .^/l -/3/(2w||) -7 2 /(16w|) (18) 
(we assume weak interactions, /3, 7 <C u>\\, hence, the oscillatios are not overdamped) and the decrement 

r = 7/2 (19) 
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of the ID monopole mode. Using the high- velocity limit Eq. ((4]) we calculate the leading terms of the two relevant 
parameters: 

/>«8^y. (20) 



v 2 Mllv) , . 



where 



(22) 



The consistency of our solution requires vN < 1 . The decrement of the monopole mode is smaller than the collision- 
induced correction to the oscillation frequency by a factor ~ v\n.(l/v). If we neglect 7, Eqs. (|18[ I2TJI 12"21 yield the 
known result [3l[ for the monopole mode frequency in a ID thermal gas. Note that dissipation has been introduced 
in [3lj phenomenologically, and up to now only the damping via essentially 3D scattering [27| has been studied. 

In the strongly interacting regime, one has to take, instead of Eq. f4|, the opposite (k r <C c) limit of Eq. (|3|). 
Then 7 decreases as c -4 as c — > 00, thus signifying integrability restoration for a longitudinally trapped gas in the 
Tonks- Girardeau limit [26j . 

We consider, as an example, a gas of AT = 10 3 87 Rb atoms trapped at T = 1 /iK in a highly elongated trap with 
the radial and longitudinal trapping frequencies uj±_ = 2n x 67 kHz and = 2ir x 10 Hz (such a strong radial 
trapping needed to satisfy the one-dimensionality conditions can be achieved in optical lattices [18]). In this case 
h 2 c/m = 2hoj±a s 32], where a s = 5.3 nm is the atomic s-wave scattering length, and we obtain i/Ri2x 1CT 4 . The 
collisional shift of the monopole mode frequency is, according to Eqs. (fT51 [2TJ)) . approximately 0.4%. The decrement 
of this mode is T s» 3 x 10~ 4 s~ 4 , which is too small to be measured experimentally. 

The method developed in the present paper applies to non-degenerate ID gases. However, by analogy with the 
results above, we can hazard a conjecture for the estimation of the decrement of the monopole mode also in the case 
of degenerate gas. Our guess is that the relation T ~ v\tv(\/v)/S.ui holds also in the latter case, but the parameter v 
has to be then estimated not from Eq. (|22|) . but instead from the the collisional shift Auj ~ vNw\\ of the monopole 
mode frequency (with respect to the frequency 2ct>ii of this mode in an ideal gas). At maximum (in the Thomas- Fermi 
regime), Auj ~ u« [33| and, hence, v ~ 1/AT. Then we can roughly estimate the decrement as V ~ hi(l/N)u)\\/N . The 
time scale on which the integrability breakdown due to longitudinal trapping can be neglected, is long compared 
to a typical duration of an ultracold-atom experiment. 



V. CONCLUSION 



We find that the combination of two effects, (i) finite time delay in wave packet propagation associated with a 
ID collision and (ii) longitudinal trapping lift the integrability of a system of identical bosons in ID with contact 
interactions. To quantify the effect of non-integrability for a non-degenerate gas, we derive a set of kinetic equations 
and apply them to the analysis of the monopole mode of a ID gas. We find this decrement to be small enough to 
be neglected in experiment and state therefore that the main source of non-integrability and thermalization in cold 
atomic Bose-gases in elongated traps in the ID regime is virtual excitation of transversal degrees of freedom that gives 
rise to effective velocity-changing 3-body elastic collisions [22T - [24l | . 

The author thanks J. Schmiedmayer and H.-C. Nagerl for helpful discussions. This work is supported by the FWF 
(project P22590-N16). 
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